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This application claims the benefit of U.S. Provisional 

Application No. 60/267,650, filed February 8, 2001 and entitled 

"In-situ Structural Health Monitoring, Diagnostics and 

p* Prognostics System Utilizing Smart Active Sensors with Embedded 

£l Intelligence and Wireless Installation," the entire disclosure 

m of which is incorporated by reference herein. 

*j£ Background Of The Invention 

l£0 It has been recognized that there is a need to monitor the 

til 

U existence and effects of damage and corrosion in structural 
U materials. Aircraft components provide a much publicized 

example, but the need exists in a variety of other structures, 
for example pipelines and building members. To detect defects, 

15 traditional modal analysis testing generally relies on 

structural excitation and vibration pickups. Structural 
excitation methods, in turn, usually employ either harmonic 
sweep or impulse excitation. The former method is more precise 
and can focus on resonant frequencies, while the latter is more 

20 expedient. Vibration pickups can measure displacement, velocity 
or acceleration. The use of self -conditioning accelerometers 
and laser velocimeters are known. Accelerometers allow 
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installation of sensor arrays that accurately and efficiently 
measure mode shapes, while lasers offer contactless measurements 
that are effective for low mass sensitive structures. 

In general, it is known to propagate electromagnetic 
5 signals through materials in order to measure anomalies in the 

materials. Ultrasonic transducers, for example, include 
, piezoelectric crystals of a certain mass. The crystals are 
S excited at a known resonant frequency so that the resonant 
f|| crystal vibration generates an ultrasonic signal. Ultrasonic 
1£R sensors, which are relatively bulky and expensive, are typically 

S5K.3! 

s used in portable devices placed in communication with the 
HI structure of interest so that the crystals propagate ultrasonic 
fl waves into the material. Sensors detecting resulting signals in 
* y the material sense patterns that indicate anomalies in the 
15 materials such as cracks or corrosion. 

Summary Of The Invention 
It is an object of the present invention to provide an in- 
situ sensor system for structural material analysis. 

This and other objects may be achieved within a method of 
2 0 detecting a damage feature in a structure. 

A plurality of thin piezoelectric ceramic sensors are 
embedded on the structure. A first sensor is excited so that 
the first sensor produces a responsive signal in the structure. 
The responsive signal is received by a second sensor, and the 
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presence or absence of at least one predetermined signal 
characteristic in the received responsive signal is determined. 
The predetermined signal characteristic is related to existence 
of the damage feature. 

The accompanying drawings, which are incorporated in and 
constitute a part of this specification, illustrate one or more 
embodiments of the invention and, together with the description, 
serve to explain the principles of the invention. 

Brief Description of the Drawings 

A full and enabling disclosure of the present invention, 
including the best mode thereof directed to one of ordinary 
skill in the art, is set forth in the present specification, 
which makes reference to the appended drawings, in which: 

Figure 1 is a cross-sectional schematic view of a thin 
piezoelectric sensor embedded on a structure for in-situ 
measurement in accordance with an embodiment of the present 
invention; 

Figure 2 is a cross-sectional diagrammatic illustration of 
a PZT wafer active sensor constrained by structural stiffness; 

Figure 3 is a graphical representation of a self -diagnostic 
performed in a sensor in accordance with an embodiment of the 
present invention; 



Figure 4 is a graphical representation of an impedance 
spectrum measured by a sensor array in accordance with an 
embodiment of the present invention; 

Figures 5a and 5b are schematic illustrations of a sensor 
array in accordance with an embodiment of the present invention 

Figure 6 is an electrical schematic view of a sensor in 
accordance with an embodiment of the present invention; 

Figures 7a and 7b are schematic illustrations of a sensor 
array in accordance with an embodiment of the present invention 

Figures 8a-8d are graphical illustrations of elastic 
measurement waves generated by a transducer array in accordance 
with an embodiment of the present invention; 

Figures 9a- 9c are schematic illustrations of embedded 
sensors in accordance with an embodiment of the present 
invention; 

Figures lOa-lOd are graphical illustrations of measurement 
and response signals resulting from use of a transducer in 
accordance with an embodiment of the present invention; 

Figure 11 is a schematic illustration of a sensor array in 
accordance with an embodiment of the present invention in 
passive detection mode; and 

Figure 12 is a block diagram illustration of a sensor and 
monitoring system in accordance with one or more embodiments of 
the present invention. 



Repeat use of reference characters in the present 
specification and drawings is intended to represent same or 
analogous features or elements of the invention. 

Detailed Description of Preferred Embodiments 

Reference will now be made in detail to presently preferred 
embodiments of the invention, one or more examples of which are 
illustrated in the accompanying drawings. Each example is 
provided by way of explanation of the invention, not limitation 
of the invention. In fact, it will be apparent to those skilled 
in the art that modifications and variations can be made in the 
present invention without departing from the scope or spirit 
thereof. For instance, features illustrated or described as 
part of one embodiment may be used on another embodiment to 
yield a still further embodiment. Thus, it is intended that the 
present invention covers such modifications and variations as 
come within the scope of the appended claims and their 
equivalents . 

Generally, the present invention relates to embedding 
relatively small transducers on a structure in order to monitor 
the structure's structural health. Several examples are 
provided, namely, active electro-mechanical impedance sensing, 
active ultrasonic sensing, and passive sensing through 
interrogation of acoustic emissions and detection of elastic 
waves indicating impacts by a foreign body. A description of an 



electromechanical (E/M) impedance damage identification strategy- 
is presented in the Appendix. Summarizing the identification 
strategy, the effect of a piezoelectric active sensor affixed to 
a structure is to apply a local strain parallel to the 
structure's surface that creates stationary elastic waves in the 
structure. The structure presents to the active sensor the 
drive-point impedance, Z 5 ^(®) = kom e (®) + c e (a)-ik e (®)/(i) . Through the 
mechanical coupling between the PZT active sensor and the host 
structure, on one hand, and through the electro-mechanical 
transductance within the PZT active sensor, on the other hand, 
the drive-point structural impedance is directly reflected into 
the effective electrical impedance as seen at the active sensor 
terminals, as shown in Figure 6. The apparent electro- 
mechanical impedance of the piezo-active sensor as coupled to 
the host structure is : 



Z(a>) = 



icoC 



/ ; ^2 Z sfr ((0) 
1 K 3] 



Z pzT (G)) + Z str (G>) y 

Where Z (co) is the equivalent electro-mechanical admittance as 
seen at the PZT active sensor terminals, C is the zero-load 
capacitance of the PZT active sensors, K 31 is the electro- 
mechanical cross coupling coefficient of the PZT active sensor 
(k 31 =d l3 /Js n s 33 ) , Z str is the impedance of the structure, and 

ZpZT IS 

the impedance of PZT active sensor. 
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The electro-mechanical impedance method is applied by- 
scanning a predetermined frequency range in the hundreds of kHz 
band and recording the complex impedance spectrum. By comparing 
the impedance spectra taken at various times during the service 
life of a structure, structural degradation and the appearance 
of incipient damage may be identified. The frequency should be 
high enough, however, that the signal wavelength is 
significantly smaller than the defect size. 

Referring to Figures 5a and 5b, consider an array of four 
active sensors. Each active sensor has its own sensing area 
characterized by a sensing radius and the corresponding sensing 
circles indicated at 18. Inside each sensing area, the sensor 
detection capability diminishes with the distance between the 
sensor and the damaged or corroded area. A damage feature that 
is placed in the sensor near field is expected to create a 
larger disturbance in the sensor response than a damage feature 
placed in the far field. Effective area coverage is insured 
when the sensing circles of several sensors overlap. The 
diagnostics of the adjacent structure is performed using the 
active (real) part of the E/M impedance. Incipient damage 
changes taking place in the structure are reflected in the 
drive-point structural impedance. The change in the structural 
drive-point impedance extensively affects the real part of the 



effective electro-mechanical impedance of the piezoelectric 
active sensor embedded on the structure. 

Figure 5a illustrates a structural crack placed in the 
sensing circle of the first active sensor. The presence of the 
crack modifies the structural field and effective drive-point 
structural impedance as seen by that sensor. At the same time, 
the crack also belongs to the sensing circle of the second 
sensor, but it is at the periphery of this circle. Thus, it is 
expected that the effective drive-point structural impedance as 
seen by the second sensor will also be affected but to a much 
lesser extent than that of the first sensor. The structural 
crack is outside the sensing circles of the third and fourth 
sensors, and their drive-point structural impedances will be 
almost unchanged. By virtue of the impedance equation above, 
changes in the drive -point structural impedance are directly 
reflected in the sensor's E/M impedance. Thus, the crack 
illustrated in Figure 5a is expected to strongly modify the E/M 
impedance of the first sensor, to slightly modify that of the 
second sensor and to leave the E/M impedances of the third and 
fourth sensors unchanged. 

Figure 5b illustrates a patch of corrosion damage placed in 
the sensing circle of the first active sensor. The corrosion 
damage also belongs, to a lesser extent, to the sensing circles 
of the second and fourth sensors. The effective drive-point 
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structural impedance seen by the first sensor will be strongly 
modified. The drive-point impedance seen by the second sensor 
will be modified to a lesser extent, and the impedance of the 
fourth sensor will be slightly modified. The drive-point 
impedance of the third sensor remains virtually unchanged. As 
indicated by the impedance equation above, these changes in 
drive-point structural impedance will be directly reflected in 
the sensors' E/M impedances. Thus, the corrosion damage 
strongly modifies the E/M impedance of the first sensor, 
somewhat modifies the E/M impedance of the second sensor, 
slightly modifies the E/M impedance of the fourth sensor and 
leaves the E/M impedance of the third sensor unchanged. 

The discussion above focuses on the use of E/M impedance 
sensing to detect defect or damage features in a structure. It 
is also possible to use active piezoelectric sensors to 
propagate elastic waves through the structure so that signal 
reflections and scattering effects may be used to detect damage 
features such as cracks and corroded areas. Although it should 
be understood that the present method may include the use of a 
single sensor, the embodiment described herein includes the use 
of a sensor array. In general, the sensors in a given array are 
clustered in a predetermined positional arrangement and are 
sequentially excited in round-robin fashion so that each sensor 
emits an ultrasonic acoustic wave into the structure. The 
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remaining sensors in the array, as well as the originating 
sensor, monitor the emitted signals and, from the received 
responses, may determine the existence and position of a damage 
feature. Excitation may be at a fixed frequency, a frequency 
5 burst or a frequency sweep. The particular frequency band may 
be selected consistently with the size of the feature to be 
identified. 

Figures 7a and 7b, illustrate an array of four active 
-it sensors in an ultrasonic sensing mode. Piezoelectric active 
1& sensors can act as both sensors and actuators. At each 

?.:s. 

i measurement, one of the four active sensors is both an actuator 

fU and a sensor, while the remaining three act only as sensors. 

O 

H;l The actuator function moves round-robin among the four sensors. 
141 Assuming the first sensor is the actuator, it generates elastic 
15 waves that propagate through the material and are sensed at the 
second, third and fourth active sensors. The properties of 
these waves are affected by the presence of impact damage or 
material defects and can be interpreted to yield damage location 
and amplitude. This method can be applied, for example, to 
2 0 detect damage arising from cracks and corrosion. 

Referring to Figure 7a, it is desirable to characterize 
crack damage in terms of its location (x c , y c ) , its size and its 
orientation (a, 9) . When a crack is present in the wave path, a 
signal generated by the actuator piezoelectric sensor will be 
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deflected, reflected and transmitted at the crack. The relative 
proportions among deflection, reflection and transmission vary 
with damage, size and orientation. In Figure 7a, the second and 
fourth active sensors receive both direct and deflected wave 
signals from the actuator (i.e. the first active sensor). The 
first sensor also acts as a receptor and detects a wave 
reflected from the crack. The third active sensor receives a 
transmitted wave having an amplitude that is a function of the 
damage size. Thus, a matrix of information in terms of event 
arrival time and amplitude may be determined from the 
information detected by the four sensors when the first sensor 
is the actuator. Thereafter, a round-robin procedure is imposed 
in which the second, third, and fourth sensors in turn become 
the actuators , An information matrix is obtained at each such 
measurement. The damage, location, size and orientation may be 
derived using conventional linear algebra or neural network 
algorithms. Generally, the use of multiple round-robin 
measurements provides more information than needed. This over- 
specification, however, mitigates data error. 

Figure 7b illustrates a corrosion damage example. 
Corrosion damage is multi-dimensional and can cover a wide area. 
In the example shown in Figure 7b, the first active sensor 
generates elastic waves that propagate through the material and 
are sensed at the second, third, and fourth sensors. The 



damaged material affects the speed and attenuation of the waves. 
As with crack damage, it is desirable to characterize corrosion 
damage in terms of location, size and orientation. Orientation 
may be defined in any suitable manner, for example the major and 
minor axes of a damage ellipse and axis inclination. In the 
example, elastic waves received by the second and fourth sensors 
do not engage the corroded area and are, therefore, unaffected. 
Waves received by the third sensor travel through the damaged 
material and are therefore modified in speed and attenuation. 
The material degradation may be so advanced that no detectable 
waves are received at the third sensor. As each of the four 
sensors emits elastic waves in round-robin fashion, an 
information matrix is stored so that the location, size and 
orientation of the damaged area are determined. Again, the 
solution may be derived from conventional linear algebra 
solutions or neural networks, and the use of multiple 
measurements mitigates errors. 

The types of elastic waves used in the damage 
identification process may vary, for example among constant 
amplitude sine waves, bursts, sweeps and impulses. Four wave- 
type examples are provided in Figures 8a - 8d. The constant 
amplitude sine wave in Figure 8a is the simplest waveform. 
Because excitation frequency should be matched with the 
structural characteristics of the damage area in order to excite 
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a structural resonance, however, use of the constant amplitude 
wave may require multiple emissions at varying frequencies to 
test for all damage types of interest. Accordingly, a frequency 
sweep as in Figure 8b may be desirable in that it permits the 
excitation of more than one frequency in the same measurement. 
An impulse signal as in Figure 8c permits the excitation of a 
wide frequency spectrum. A frequency burst as in Figure 8d may 
be desirable in that it includes a dominant frequency that can 
be tuned to particular structural and flaw-size requirements. 
Since its frequency content is known, a burst wave can be 
readily detected and filtered from background noise, and its 
limited duration facilitates the identification and analysis of 
burst reflection from defects and boundaries. 
Experimental Results 

Lamb-wave propagation experiments were conducted on aging 
aircraft panels using a number of PZT active sensors affixed at 
various locations. Experiments were performed to verify the 
wave propagation properties and to identify reflections due to 
intrinsic construction features of the panels, for example, 
rivets, splice joints, etc. The illustrations in Figures 9a - 
9c provide a crack detection example. Referring to Figure 9a, a 
series of rivets 20 attach to panel sections 22 and 24 so that a 
side edge 23 of panel 24 faces toward the right in the view 
shown in Figure 9a. A piezoelectric active sensor 26 is 
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attached by adhesive to panel 24 at a distance of 200mm from 
rivets 20. A more complicated array of rivets is provided on 
panel 24 in Figure 9b. The panel shown in Figure 9c includes 
the same rivet arrangement, with the addition of a simulated 
crack (12.75mm EDM hairline slit) 28 starting from the right- 
most rivet at the top horizontal row. That is, Figure 9c 
illustrates the damaged panel. Experiments were also conducted 
on a panel (not shown) having no rivets or cracks. 

In the clear panel (not shown) and the panels of Figures 9b 
and 9c, PZT sensors 2 6 emitted an approximately 2 5mV continuous 
frequency wave. Figure 10a shows the emitted and received 
signals on the undamaged panel without rivets (not shown) . The 
graph shows the emitted signal, centered at around 5.3 micro- 
seconds, and multiple reflections from the panel edges. The 
reflections start to arrive at approximately 60 micro- seconds . 
Figure 10b illustrates the signals for the panel shown in Figure 
9b. The signal features reflection from the rivets and multiple 
reflections from the panel edges. The reflection from the 
rivets arrives at approximately 42 micro- seconds , indicating an 
approximate time-of -flight having the value TOF = 37 micro- 
seconds. This TOF is consistent with traveling 200mm at a group 
velocity of approximately 5.4km/s. Figure 10c shows the signal 
recorded on the damaged panel of Figure 9c. The signal response 
is similar to that shown in Figure 10b, but is stronger at the 
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42 micro-second position. Subtracting the signal of Figure 10b 
from that of Figure 10c, as shown in Figure lOd, identifies the 
effect of the crack. That is, the strong wave pack centered on 
42 micro-seconds illustrates the response from the crack. 

Accordingly, where an array of sensors is disposed on or in 
a panel, and where signals are initially measured in the 
undamaged structure to obtain baseline measurements, the 
baseline measurements may be subtracted from subsequent 
measurements to identify the existence of damage features. The 
time -of -flight determines the distance from the respective 
sensor. By measuring distances from multiple sensors in the 
array, the location of a damage feature with respect to the 
sensors may be determined. Moreover, the use of multiple 
measurements, sequentially activating the individual sensors in 
the array in round-robin fashion, provide for precise location. 

The examples provided above, i.e., impedance sensing and 
ultrasonic sensing, rely on active sensors to measure structural 
characteristics. The present invention may also, however, be 
used in passive modes to detect structural damage. In certain 
of there embodiments, a plurality of sensors is disposed in a 
predetermined orientation relative to each other and at known 
positions on the structure. The sensor outputs are monitored 
intermittently or continuously, even though the sensors may not 
be engaged in either of the active measurement procedures 



described above. Damage events may be identified through the 
reception of stress waves generated in the structure through 
impacts or other material disruptions. Certain waves may, for 
example, indicate an occurrence of a low-velocity impact. The 

5 sensor may also, however, detect acoustic emission signals that 
indicate damage has occurred. By determining and recording the 
location and time of damage events, a record may be compiled to 

S predict the structure's remaining operative life. 

fair 
i! 

f|| Figure 11 illustrates a damage event taking place in a zone 

lip; monitored by four piezoelectric transducers acting as passive 
s sensors. In this example, the damage under consideration is a 
[M low-velocity impact. The impact location and intensity are 
determined from elastic wave signals received by the four 
sensors. The arrival time at each sensor is proportional to the 
15 distance between the sensor and the impact location. 

Accordingly, a simple triangulation is used to determine the 
impact's position with respect to the sensors. This information 
is provided to a computer (not shown) that time-stamps the 
information so that a damage event record is kept. With four 
2 0 sensors, more data is available than necessary. Thus, an 

optimal solution algorithm, for example least squares, is used. 
Damage intensity is determined simultaneously by inclusion in 
the variable list. A neural network algorithm may be applied 
instead of explicit calculation. 
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Information-rich data received from the sensors is 
preferably processed such that only relevant content is 
retained. The elimination of spurious noise and far- field 
disturbances, and the identification of useful information 
5 directly related to structural damage, is addressed through 
appropriate signal processing methodologies. Mathematically, 
the determination of the structure's physical condition using 
pi sensor measurements is a nonlinear inverse problem. Signal 
ffj processing and interpretation methods may be application- 
3JI specific (model -based) or generic (non-model) . Conventional 

itJS 

p (for example, Fourier analysis) and advanced (for example, 

III Wavelet analysis and digital filters) algorithms available in 

P 

y| specialized software packages (for example, MATLAB, LABVIEW, 

MATHCAD, or AUTO S I GNAL ) may be used, as well as neural networks 

15 and expert systems. 

In operation, transducers in a measurement array execute 
self -diagnostics (as described in the Appendix) prior to damage 
detection measurement cycles. Temperature and humidity are the 
primary environmental changes that affect sensor calibration. 

2 0 Generally, these changes may affect the sensor itself and/or its 
adhesive interface with the structure. Although the values of 
the piezoelectric coefficients vary with temperature, the 
general aspect of the frequency response curve is only shifted. 
Thus, by calibrating the frequency shift over a variety of 
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temperatures, the system may compensate for temperature changes, 
thereby maintaining the general nature of the E/M impedance 
spectrum. The sensor/structure adhesion may be affected by both 
humidity and temperature, and the combination of humidity and 
temperature change may be especially detrimental to the adhesive 
interface. Thus, the self -diagnostic described in the Appendix 
with respect to Figure 3 may be employed to identify disbonded 
sensors . 

Each sensor may be formed as a small square or disk of thin 
piezoelectric ceramic material. In one embodiment, the length 
and width of each wafer is less than 13mm, +/- 0.13mm, and 
between 0.2mm and 0.49mm thick, +/- 0.025mm, respectively. 
Electrodes are attached to the sensors, which communicate with a 
personal computer, other computing devices and/or other system 
components, through hard-wire or wireless connections. Sensors 
may be embedded on structures by, for example, insertion into a 
composite material during the layering process so that the 
sensor is actually within the structure or by adhesive 
attachment to the structure's outer surface as discussed above. 

Figure 12 provides an overview of an embodiment of the 
present invention. System 3 0 includes one or more structures 32 
having one or more transducer arrays disposed thereon as 
described above. For ultrasonic acoustic wave sensing, an 
ultra-high frequency high-energy tone-burst source 34 excites 
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individual sensors within each transducer array in round-robin 
fashion through a programmable wide-band bi-directional twenty- 
channel scanner/selector 36. These components communicate with 
a personal computer 3 8 through an ultra- fast data acquisition 
module 40. Computer 38 also acquires measurement data from the 
transducer array through data acquisition module 40. 

In wireless applications, the system may include a wireless 
communication data-burst transmitter/receiver server 42 that 
excites and reads the transducer array through a plurality of 
data-concentrators 44. That is, assuming there are a plurality 
of transducer arrays disposed about a structure, each transducer 
is attached to an antenna that communicates with a data- 
concentrator 44 dedicated to the transducer's array. The one or 
more data-concentrators communicate with transmitter/receiver 
server 42 which, in turn, wirelessly communicates with an 
antenna associated with computer 38 or acquisition module 40. 

For impedance sensing, an adjustable-voltage-power gain- 
phase impedance analyzer 46 excites the transducers, which 
output measurement information back to analyzer 46. Analyzer 46 
includes software algarithons to analyze the sensor data to 
determine the locations and orientation of damage features as 
described above. Alternatively, analyzer 46 may forward the 
data to computer 38 for analysis. Generally, computer 38 houses 
the system's software components, the operation of which is 
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generally described above and which may include non-destructive 
evaluation and imaging software package 48, monitoring 
diagnostics and analysis software package 52 and/or artificial 
intelligence, neural -network and data mining software 54. 

While one or more preferred embodiments of the invention 
have been described above, it should be understood that any and 
all equivalent realizations of the present invention are 
included within the scope and spirit thereof. The embodiments 
depicted are presented by example only and are not intended as 
limitations upon the present invention. Thus, it should be 
understood by those of ordinary skill in this art that the 
present invention is not limited to these embodiments since 
modifications can be made. Therefore, it is contemplated that 
any and all such embodiments are included in the present 
invention as may fall within the literal or equivalent scope of 
the appended claims. 
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Because of their intrinsic electro-mechanical (E/M) 
coupling, commercially available low-cost piezoceramics can be 
used as both sensors and actuators . The frequency bandwidth of 
these devices is orders of magnitude larger than that of 
conventional shakers and impact hammers. Small piezoelectric 
ceramic wafers can be permanently attached to structural 
surfaces to form sensor and actuator arrays that permit 
effective modal identification in a wide frequency band. 
Assuming a single degree of freedom, the electrical admittance, 
as measured at the terminals of the piezoelectric (PZT) wafer 
attached to a structure, is as follows: 



where C is the electrical capacitance of the PZT active sensor, 
Z s tr((£>) is the one-dimensional structural impedance as seen by 
the sensor, and Z A (co) is the quasi-static impedance of the 
sensor. At coupled- system resonance, the real part of the E/M 
admittance has a distinct peak. However, due to the additional 
stiffness contributed by the PZT wafer, the system's natural 
frequency may shift, although this effect diminishes greatly as 
the difference in size between the wafer and the structure with 
which it is used increases. 




(1) 



The E/M impedance/admittance response may be used to 
measure the structural dynamics of a material . Referring to 
Figure 1, consider a one-dimensional structure (for example a 
beam) 10 upon which is attached a ceramic PZT active sensor 12 
having a length l a , at a position on beam 10 between x a and x a + 
2 a . Upon activation, the PZT active sensor expands by epzt. This 
generates a reaction force F PZT from the beam onto sensor 12 and 
an equal and opposite force from the sensor onto the beam. This 
force excites the beam. At the neutral axis, the effect is felt 
as an axial force excitation, N PZT/ and a bending moment 
excitation, M PZT . As the active sensor is electrically excited 
with a high-frequency harmonic signal, it induces elastic waves 
into the beam structure. The elastic waves travel sideways into 
the beam structure, setting the beam into oscillation. In a 
steady-state regime, the structure oscillates at the PZT 
excitation frequency. The reaction force per unit displacement 
(dynamic stiffness) presented by the structure to the PZT 
depends on the internal state of the structure, on the 
excitation frequency, and on the boundary conditions: 

where u PZT {®) is the displacement amplitude at frequency co, F P2T (co) 

is the reaction force, and Jc str (co) is the dynamic stiffness. The 
symbol A signifies amplitude. Since the size of the PZT is very 
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small with respect to the size of the structure, equation (2) 
represents a point-wise structural stiffness. 

Consider a PZT wafer of length l a/ thickness t at and width 
b a , undergoing longitudinal expansion, u lr induced by the 
thickness polarization electric field, E 3 . The electric field is 
produced by the application of a harmonic voltage V(t) = Fe icot 
between the top and bottom surfaces (electrodes) . The resulting 
electric field, E = V/t, is assumed uniform with Xi {dE/dx^o) . 
The length, width, and thickness are assumed to have widely 
separated values ( t a <<i> a <<l a ) such that the length, width, and 
thickness motions are practically uncoupled. 

Relevant equations regarding the PZT material are: 

S^sfc+d^ (3) 

A = <W1 +4^3 ( 4 > 

where S ± is the strain, T ± is the stress, D 3 is the electrical 
displacement (charge per unit area) , s* is the mechanical 
compliance at zero field, e$ is the dielectric constant at zero 
stress, d 3 i is the induced strain coefficient, i.e., mechanical 
strain per unit electric field. 

Since the piezoelectric active sensor is ultimately bonded 
to the structure, an elastically constrained sensor is 
considered. As a result of this assumption, the active sensor 



is constrained by structural stiffness, k stri as shown in Fig. 2 
For symmetry, the structural stiffness has been split into two 
end components, each of size 2k str . The boundary conditions 
applied at both ends connect the resultant of internal stresses 
with the spring reaction force, i.e., 

2l(H)V.=-2**r«l(H) (5) 
T^\l a )b a t a =2k str u\-\l a ) 

Substitution of Eq. (5) into Eq. (3) gives: 

«i'(H)v fl =-2^ r 4«,(H)+^3 

? (6) 

"i'(-H)V«=2^^-« 1 (-i/ a )+^A 

Introducing the quasi -static stiffness of the PZT active sensor, 



and the stiffness ratio 



h - A « (1) 

K PZT - Ei ^ 1 ' 



r = ^ (8) 

"■PZT 



Now Eq. (6) can be re -arranged in the form: 

< (±H)*A ±T7-"l [¥a] = d i\ E i ( 9 ) 



Using Newton's law of motion, r/ = p«, , and the strain- 
displacement relation, Eq. (3) yields the axial waves 
equation: 

«,=c 0 V (10) 
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where Q=d()/dt , and ()'=£()/& , while c 2 a = \/ Q s* is the sound speed. The 
general solution of Eq. (10) is: 

u 1 (x,t) = u 1 (x)e' s>t where u ] (x) = (C l sinyx + C 2 cosyx) (11) 

The variable y = co/c« is the wave number, and O signifies the 
harmonic motion amplitude. The constants C ± and C 2 are to be 
determined from the boundary conditions. 

Substitution of the general solution (Eq. 11) yields the 
following linear system in d and C 2 : 



((p coscp + r sin<p ) C x - (9 sincp - r coscp ) C 2 = j u l 
(<pcos<p + rsm(p)C { + (<psin(p - rcos<p)C 2 = ^ 



(12) 



where u 



l isA=^A' l a an d 9=iY^- Upon solution, 




smyx 



(13) 



Equation (4) can be re-expressed as: 




(14) 



where 



= 4/(^3) is the E/M coupling factor. Integration of Eq. 



(14) yields the charge: 



1 2 t [ L 



2 2 




J 



(15) 



where z^bjjt^c is the conventional capacitance of the 



piezoelectric wafer. For harmonic motion, i=m-Q . Recalling the 



expressions Y = l/v and Z=T X for the electric admittance and 
impedance, the admittance and impedance expressions for a PZT 
active sensor constrained by the structural substrate with an 
equivalent stiffness ratio r are: 



y = i©-c 



1-k 



1 



(pcot<p + r 



Z = 



\g>C 



1-- 



<pcot<p + r J 



(16) 



In this equation, the structural stiffness ratio, r, is 
additive to the sensor resonance term, (pcotq> . When the PZT 
active sensor is used in a frequency sweep, the apparent 
structural stiffness, k str , will vary with frequency, going 
through zero at structural resonances and through extreme values 
at structural anti -resonances . Eqs. (16) infer that both 
structural resonances and sensor resonances will be reflected in 
the admittance and impedance frequency spectra. 

Analysis of the asymptotic behavior of Eqs. (16) reveals 
important facts. , As r 0, i.e., vibrations of a free- free 



sensor, 



(pco\(p 



As r becomes infinitely large, 



i.e., vibrations of a clamped sensor, Ydamped m ' c [ l 



On the 



other hand, as yla 0 (i.e., quasi-static sensor conditions) 



The expressions contained in Eqs . (16) bridge the gap 
between high-frequency sensor- focused analysis and low- frequency 
structure-focused analysis. The present results cover the 
complete frequency spectrum and encompass both structure and 
sensor dynamics. 

The effect of structural and sensor damping can be easily 
introduced in Eqs. (16) by the use of complex notations: 



where r is the frequency-dependent complex stiffness ratio that 
reflects the structural point -wise dynamics and the sensor 
dissipation mechanisms. 

The response of the structural substrate to the PZT 
excitation is deduced from the general theory of beam 
vibrations. The excitation forces and moments acting upon the 
beam structure are derived from the PZT force, F PZT =F PZT e im , using 
the beam cross-section geometry: 



The space-wise distribution of excitation bending moment 
and axial force are expressed using the Heaviside function, H(x 
~ x a) , defined as H(x - x a ) = 0 for x < x a , and H{x - x a ) -1 for 





(18) 



Ma=Fpzr\, N a — F PZT (19) 



KixJ) = N a [-H{x-x a ) + H{x-x a -l a )\e m (20) 
K(x,t) = -M a [-H(x-x a ) + H(x-x a -l a )}e"»< (21) 
Equations (2 0) and (21) correspond to axial and flexural 
vibrations, respectively. Axial vibrations modes are usually of 
much larger frequency than flexural vibration modes, but their 
vibration frequencies are commensurable with those of the PZT 
active sensors. Thus, both axial and flexural vibrations are 
considered. 

The equation of motion for axial vibrations is: 

pA ■ u(x,t) - EA • u\x,t) = N f e (x,t) (22) 
Assume modal expansion 

= (23) 

where X n (x) are orthonormal mode shapes, i.e., \x m x n dx = § mn , with 
8^=1 for m = n, and 0 otherwise. C n are the modal amplitudes. 
Substitution of Eq. (20) into (22) yields: 

pA ■ ii(x,t) - EA ■ u*(x,t) = N a [-S(x-x a ) + S(x-x a -l a )]-e mt (24) 

where 5 is Dirac's function. Since the mode shapes satisfy the 
free-vibration differential equation EA-X" n +($* *pA-X n = 0 , 
multiplication by X n (x) and integration over the length of the 
beam yields : 

c ^^h-^- x ^ +x ^ +l ^ (25) 
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Hence , 

«(*,t) = ^± ~ X ' {x ' ) l X '{ X ' +I ' ) -X,(x). t t'> (26) 
pA „ =0 co n -co 

For Euler-Bernoulli beams, the equation of motion under 
moment excitation is: 
5 pA- w(xj) - EI • w"\xj) - ~M n e {x,t) (27) 

Substitution of Eq. (21) into (27) yields: 

:; : 
z'z 

H pAw(x,t) + EIw""(x,t) = M a [-S'(x-x a ) + S'(x-x a ~l a )]e m ' (28) 

r$ where 8' is the first derivative of Dirac's function (5' = //"). 

ssisis 

tt Assume modal expansion w(x,t)= £ C„X n (x)-e"" , where X n (x) are the 

i y »=/s 

£| orthonormal bending mode shapes, and N lf N 2 are the mode numbers 

£9 

ftl enduring the frequency band of interest. Since the mode shapes 
satisfy the free-vibration differential equation El-x;"=«> 2 n - P A-x n , 
15 multiplication by X n (x) and integration over the length of the 
beam yields: 

C ^^-^l X ^'( x -^- b '( x - x ^0¥ X (29) 

Integration by parts and substitution into the modal 
expansion expression yields: 

Hence, 
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J xt)= M a j, -X'„(x a ) + X' a (x a+ l a ) 
pA% G>l-a> 2 



(31) 



15 



To obtain the dynamic structural stiffness, k str , presented 
by the structure to the PZT, the elongation between the two 
points, A and B, connected to the PZT ends is first calculated. 
Simple kinematics gives the horizontal displacement of a generic 
point P placed on the surface of the beam: 

u p (t) = u(x)~W(x) , (32) 

where u and w are the axial and bending displacements of the 
neutral axis. Letting P be A and B, and taking the difference, 
yields : 

»iw(0 = «a(/)-»4<0 = «(*.,0^^ (33) 

Using Eqs . (19), (26), and (31), the amplitude of Eq. (33) 
becomes 



pA 



x h^^ [w' n {x a+ l a )-W' K {xS$ 
2jf ' 



(34) 



where differentiation between axial and flexural vibrations 
frequencies and mode shapes was achieved by the use of n u , ©„ , 

U »A X ) and n «" ' w k( x ) > respectively. Dividing Eq. (34) by F PZT 
yields the structural frequency response function (FRF) to the 
20 Single Input Single Output (SISO) excitation applied by the PZT 
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active sensor. This situation is similar to conventional modal 



testing with the proviso that the PZT wafers are unobtrusive and 
permanently attached to the structure. For convenience, the 
axial and flexural components of the structural FRF are 
expressed separately, i.e., 

„.(„>. J^fc^^kM, ,35, 



^^[^[WKM (36) 



Modal damping, was introduced to provide practical 

veridicality to the model. The FRF' s are additive, and the total 
FRF is simply 

//(©) = H u (g>) + H w (g>) (37) 

The SISO FRF is the same as the dynamic structural 

compliance, as seen by the PZT wafer active sensor. The dynamic 
structural stiffness is the inverse of the structural 
compliance, i.e., 



U PZT 



[tU*.+/ a )-M*.)T 



(38) 



For free-free beams axial and flexural components [31] : 
U„ t (x) = A„cos(y K x), A^-m , Y % =?f- , <*„,=y % c , c = ^ , n u = 1,2,... (39) 



Numerical values of /-y„ and a % for n w <5 can be found in 
Blevins, R.D., (1979) "Formulas for Natural Frequency and Mode 
Shape," page 108, Litton Educational Publishing Inc.; for 5<n w , 

r„ w =j 2 ' and a n ^\ . 

5 Experimental Results 

The analytical model described above was used to perform 
p several numerical simulations that directly predict the E/M 
%l impedance and admittance signature at an active sensor's 
P terminals during structural identification. Subsequently, 
experiments were performed to verify these predictions. The 

e 

J:| simulation conditions identically represent specimens consisting 

Q 

|| of small steel beams (E = 200GPa, p = 7750 kg/m 3 ) of various 

o 

l|| thicknesses and widths. All beams were 100mm long with various 
widths, bi = 8mm (narrow beams) and b 2 = 19.6mm (wide beams) . 

15 The nominal thickness of the specimen was hi = 2.6mm. Double 
thickness specimens, h 2 = 5.2mm, were created by gluing two 
specimens back- to-back . Thus, four beam types were used: 
narrow-thin, narrow- thick, wide-thin, and wide-thick. The 
comparison of wide and narrow beams was aimed at identifying the 

2 0 width effects in the frequency spectrum, while the change from 
double to single thickness was aimed at simulating the effect of 
corrosion (for traditional structures) and 

disbonding/delamination of adhesively bonded and composite 
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structures. All specimens were instrumented with thin 7mm 
square PZT active sensors (1 = 7mm, b = 7mm, and t = 0.22mm) 
placed at x a = 4 0mm from one end. 

The numerical simulation was performed using a MATHCAD- 
coded simulation based on modal analysis theory, assuming a 
damping coefficient of 1% for steel beams. Numerically exact 
expressions for the axial and flexural frequencies and mode 
shapes were used. The simulation was performed over a modal 
subspace that incorporates all modal frequencies in the 
frequency bandwidth of interest. The theoretical analysis 
indicates that these frequencies should be identical with the 
basic beam resonances. The "calc" columns of Table 1 show the 
first six predicted resonances for axial and flexural 
vibrations. During the experiments, a Hewlett-Packard 4194A 
Impedance Analyzer was used to record the E/M impedance real 
part spectrum in the l-3 0kHz range. To approximate the free- 
free boundary conditions, the beams were supported on common 
packing foam. The beam natural frequencies were identified from 
the E/M impedance spectrum. The results are given in the U EXP" 
columns of Table 1. When the beam thickness was doubled, the 
frequencies also doubled. This is consistent with theoretical 
prediction. The error between theory and experiment, however, 
seems larger for the double thickness beam. 



The experimental results indicate frequency clusters, which 
move to higher frequencies as beam width is reduced. The 
frequency clusters appear to be associated with width 
vibrations. Width vibrations are also influenced by thickness, 

5 i.e., they shift to higher frequencies as thickness increases. 
This is also noticeable in Table 1, which shows that the lowest 

u cluster appeared for the single- thickness wide specimen and that 

P| the highest cluster appeared for double- thickness narrow beam. 

i 

flj For purposes of comparison, dynamic identification was also 

l4* attempted with conventional modal analysis methods. A small 
^ steel beam with dimensions identical to the "narrow thin" 
^ specimen in Table 1 was instrumented with two CEA-13 -240UZ-12 0 
~? strain gauges connected in half bridge configuration to a P-3500 
strain indicator available from Measurements Group, Inc. The 
15 specimen was suspended in a free- free configuration and excited 
with a sharp impact. The resulting signal was collected with an 
HP 54601B digital oscilloscope and numerically processed on a 
personal computer. Standard signal analysis algorithms (FFT) 
were used to extract the frequency spectrum. The first natural 
20 frequency (1.387kHz) was clearly displayed. The second natural 
frequency (3.789kHz) could also be identified but with a much 
weaker amplitude. These results were consistent with 
theoretical values and the experimental results presented in 
Table 1, but the impact excitation method was not able to excite 



34 



other higher frequencies depicted in Table 1, most probably due 
to bandwidth limitations. 

In order to obtain consistent experimental results, a set 
of sensors was subject to a calibration procedure consisting of 
geometrical and electrical measurements and measurements of the 
intrinsic E/M impedance and admittance spectra of the PZT active 
sensors. The results show small dimension variations from 
sensor to sensor. The mean and standard deviation values for 
length/width and thickness were 6.95mm, +/- 0.5%, and 0.224mm, 
+/- 1.4%, respectively. Electrical capacitance was measured on 
an in-process quality check to be applied during each step of 
sensor development and during the sensor installation process. 
Mean and standard deviation values were 3.276mF, + /- 3.8%. 

The intrinsic E/M impedance and admittance spectra of the 
PZT active sensors, before being attached to the structure, were 
measured with an HP 4194A Impedance Phase-Gain Analyzer. The 
PZT wafer was centered on a bold head and held in place with a 
probe tip. Thus, center clamping conditions were simulated, and 
the wafer could perform free vibrations while being tested. 
Mean values of 2 51kHz and standard deviations of +/- .2% were 
obtained. 



Table 1 



Theoretical and experimental results for wide and 
narrow beams with single and double thickness 



# 


Beam #1 (narrow thin) 


Beam #2 (Narrow thick) 


Beam #3 (wide thin) 


Beam #4 (wide thick) 


Calc. 
KHz 


Exp. 
kHz 


A% 


Calc. 
kHz 


Exp. 
kHz 


A% 


Calc . 
kHz 


Exp. 
KHz 


A% 


Calc. 
kHz 


Exp. 
kHz 


A% 


1 


1.396 


1.390 


-0.4 


2.847 


2.812 


-1.2 


1.390 


1.363 


-1.9 


2.790 


2.777 


-0.5 


2 


3.850 


3 .795 


-1.4 


7. 847 


7.453 


-5.2 


3.831 


3.755 


-2 


7.689 


7.435 


-3.4 


3 


7.547 


7.4025 


-2 


15.383 


13 .905 


-10.6 


7.510 


7.380 


-1.7 


15.074 


13 . 925 


-8.2 


4 


12.475 


12.140 


-2.7 




29.650 




12.414 


12.093 


-2.6 




21.825 




5 


18 . 635 


17.980 


-3.6 


25 .430 


21 . 787 


-16 . 7 


18.545 


17 . 965 


-3.2 


24 .918 


22 .163 


-12.4 


6 




24 .840 












24 . 852 










7 


26.035 


26.317 


1 


26 .035 


26.157 


0.5 


26.022 


26.085 


0.2 


25 .944 


26.100 


0.6 






Cluster 
175 kHz 






Cluster 
210 kHz 






Cluster 
35 kHz 






Cluster 
60 kHz 





□ The active sensors used in the experiments were very small 

f|| and therefore did not significantly disturb the dynamic 
4* properties of the structure under consideration. Table 2 

;L presents the mass and stiffness for the sensor and structure. 

f 1 1 

ljftj For purposes of comparison, the table also presents the mass of 

1% an accelerometer. 

m 

Table 2 

15 Numerical illustration of the non- invasive 

properties of the piezoelectric wafer active sensors 



Sensor type 


Mas s , g 


% of 
structural 
mass 


Stiffness, 
MN/m 


% of 
structural 
stiffness 


Active sensor 


0.082 


0.5% 


15 


1.5% 


Structure 


16.4 


N/A 


1000 


N/A 


Accelerometer 
:352A10, PCB 
Piezot ronics 


0.7 


4.3% 


N/A 


N/A 



Table 2 illustrates the non-invasive properties of 



20 piezoelectric wafer active sensors. The mass and stiffness 

editions of the sensors are within about 1% (0.5% for mass, and 
1.5% for stiffness). 
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As described in more detail above, the present invention 
contemplates embedding piezoelectric wafer transducers on 
mechanical structures for potentially long periods of time. 
Thus, it is desirable to intermittently monitor the integrity of 
a transducer or transducer array. Accordingly, a self- 
diagnostic may be employed in which the reactive part of the 
transducer's impedance is measured over a relatively broad 
frequency range. Piezoelectric active sensors are predominately 
capacitive devices dominated by reactive impedance. Baseline 
impedance signatures can be used to identify defective active 
sensors. Referring to Figure 3, for example, a reactive 
impedance spectrum 14 of a well -bonded PZT sensor is compared 
with an impedance spectrum 16 of a disbonded (free) sensor. The 
appearance of sensor free -vibration resonance and the 
disappearance of structural resonances provide un-ambiguous 
features that indicate the degree of the transducer's integrity. 

Piezoelectric active sensors, and the associated structural 
dynamic identification method described above based on electro- 
mechanical impedance response, is well -suited for small rigid 
machinery parts having natural frequencies in the kHz range. 
For example, Figure 4 illustrates the natural frequencies 
identified by two PZT active sensors disposed on an aircraft 
turbo-engine blade, one on the blade's main portion and one at 
the blade root. Following identification of such impedance 



spectra from the sensors, the sensors may thereafter be 
periodically monitored for deviations in the spectra that could 
indicate structural defects in the blade. 
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